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Abstract. Liquid flow propelled by capillary forces is one of the most important 
-^.j transport mechanisms in porous environments. It is governed by a fascinating interplay 

Q of interfacial, viscous drag as well as gravitational forces which liquids encounter upon 

'^ invasion into geometries with often complex topologies, such as capillary networks of 

"t^ trees or interconnected fractures in soils and ice. Whereas a detailed understanding of 

C^ this phenomenon has been achieved on macroscopic scales, its phenomenology is poorly 

I explored on the meso- and nanoscale, where predictions regarding its dynamics are 

' ^ hampered due to a possible break-down of continuum hydrodynamics and conflicting 

pi. reports with respect to the conservation of the fluidity and capillarity of liquids upon 

(^ spatial conflnement. Here, we present fundamentals, concepts and an experimental, 

gravimetric study on the capillarity-driven invasion dynamics of liquids in networks of 
T-H pores a few nanometers across in monolithic, nanoporous silica glass (porous Vycor). 

^ A variation of the complexity of the building blocks of the liquids investigated (water, 

~r chain- like n-alkanes, silicon oils and rod-like liquid crystals) along with a variation 

r- — of the humidity and the temperature upon spontaneous imbibition allows us to 

^^ gain information regarding the fluidity and capillarity of liquids in such nanoporous 

l/~j environments. We observe square-root of time imbibition dynamics for all liquids 

^P applied, which we can quantitatively describe by both a conserved bulk fluidity in 

the pore center and bulk capillarity at the advancing menisci, if we assume a sticky 

boundary layer (negative velocity slip length). Moreover, pecularities of nanopore- 

confined liquids, such as transport via the vapor phase leading to preadsorbed liquid 

j^ layers, have to be properly accounted for. Upon increasing the chain-length in the case 

of the n-alkanes, we found hints towards a transition from stick- to slip-flow at the pore 
walls with increasing chain-length and thus polymeric behavior. Meniscus freezing, 
that is the formation of a surfactant-like, rectified monolayer at the advancing menisci 
of the imbibition front, is reported for n-tetracosane confined in porous Vycor. For 
the rheology of a rod-like liquid nematogen (80CB) we found no hints of the viscosity 
drop upon entering into the nematic phase, typical of the bulk rheology of this liquid 
crystal. 



PACS numbers: 47.55.nb, 47.61.-k, 87.19.rh, 47.57.Lj, 47.56.-Fr 
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1. Introduction 

Spontaneous imbibition, that is the capillarity-driven invasion of a liquid in porous 
material is a very common phenomenon. From everyday life almost everybody knows 
the effect of coffee or tea being absorbed by a cube of sugar. A rather unpleasant 
instance of being confronted with its consequences is rising moisture in your basement 
walls. Another well-known example is the transport of water from the roots of a tree 
up to its limbs, branches and leaves, which is partly driven by capillary action [T]. 

There has been significant progress with regard to a quantitative understanding of 
this phenomenon, most prominently with regard to the transport through rocks and 
soils, starting with the seminal work of the French engineer Henry Darcy (1803-1858) 
[21 E]. Almost all experimental and theoretical studies, published to date, deal however 
with porous structures where the characteristic length scales, meaning the typical pore 
diameters, are on the macroscale, that is much larger than a typical size of the basic 
building blocks of the flowing liquids [1] . 

After a short review of, and introduction into, the imbibition phenomenology of 
liquids in porous media in general, we will present a study on the spontaneous imbibition 
of liquids in pores a few nanometers across. We will vary the complexity of the building 
blocks of the liquids from water to simple short-chain n-alkanes and silicon oils to 
rod-like liquid crystals. Examinations on the phase transition behavior of nanopore- 
confined liquids will be part of this study as well. The goal is to scrutinize which 
principles regarding the rheology of these liquids, known from the bulk state, can be 
safely transferred to the transport across this nano- or mesoscopic geometries. 

Given the emerging interest in micro- and nanofiuidic applications [5j such a 
study on the non-equilibrium behavior of liquids in extreme spatial confinement, 
on the nanometer scale, is not only of fundamental interest, but also of rather 
practical importance [B] . Moreover, the increasing use of nanoporous matrices as hard 
templates |7j for the preparation of well-defined nanoscopic, soft-matter structures, 
such as nanotubes and nanorods |H], necessitates a profound understanding of this 
phenomenology. 

2. Fundamentals of liquid imbibition in nanoporous solids 

2.1. Liquid flow in isotropic pore networks: Darcy 's law 

In order to introduce the fundamentals of liquid flow in nanoporous media, we start 
with a much simpler phenomenology: The flow of a liquid through a pipe is a prime 
example for the direct application of the Navier-Stokes equation resulting in the famous 
law of Hagen- Poise uille. For a given pressure difference Ap applied along a cylindrical 
duct with radius r and length i the volume flow rate V is determined by 

4 

vr r 
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Here r) denotes the dynamic viscosity of the flowing hquid. In the following we will 
present a simple concept that allows one to extend the applicability of Eq. ([I]) to the 
flow through a complex pore network. 

2.1.1. Characterization of a pore network Porous Vycor glass (code 7930) provided by 
Corning Incorporated was used over the course of this study. It is produced through 
metastable phase separation in an alkali-borosilicate system followed by an extraction 
of the alkali-rich phase. The final glass consists of a sponge-like network of tortuous 
and interconnected pores with mean pore radii on the nanometer scale embedded in a 
matrix mostly consisting of Si02 [H] • 

In general such an isotropic pore network can be characterized by three quantities. 
The mean pore radius tq and the volume porosity 0o are probably the most intuitive ones 
among them. In the course of this study two different batches of the monolithic glass 
were applied that differed in tq whereas they coincided in 0o ~ 0.3. For convenience we 
will refer to them as V5 {tq = 3.4 nm) and VIO (rg = 4.9 nm) from now on. These matrix 
properties were accurately ascertained employing nitrogen sorption isotherms conducted 
at 77 K and through a subsequent analysis within a mean field model for capillary 
condensation in nano- and mesopores. More details on the matrix characterization and 
preparation are given in Ref. [TD]. 

With only these two parameters a porous cuboid with edge length a (and cross- 
sectional area A = a^) consisting of 



,A / Koid nvrr^a 

-^ ^ 00 = T-r = 2 

TT Tn \ ^sample Q' 



^ = ^^ ^ 00 = T} = o (2) 



cylindrical pores with radius tq and length a can be constructed. Assuming the 
capillaries to be aligned in flow direction the flow rate through the whole matrix is 
then given by n times the single pore flow rate Eq. ([I| with r = tq and i = a. However, 
so far this description still lacks information on the orientation of the pores. 

To account for the isotropy of the network as indicated in Fig. Q^(left) it is necessary 
to introduce a third parameter, the so-called tortuosity r along with the transformation 

V ^ -V (3) 

T 

of the volume flow rate Eq. (fTl). Pores totally aligned in flow direction would yield t = 1, 
whereas isotropic distributed pores would result in r = 3. This can be seen very easily. 
For a random orientation only every third pore is subjected to the pressure gradient 
and hence contributes to the flow. Therefore the net flow rate has to be divided by the 
factor three. 

To date serveral techniques have been applied to extract the tortuosity of the 
isotropic pore network in Vycor glass. Deducing the diffusion coefficient of hexane 
and decane by means of small angle neutron scattering (SANS) measurements r was 
found to be in the range of 3.4 - 4.2 [TT]. Gas permeation measurements performed 
with an in-house apparatus resulted in r = 3.9 ± 0.4 p^. Finally, calculations based on 
three-dimensional geometrical models yielded a value of approximately 3.5 
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Interestingly, all values show a significant deviation from r = 3 as derived from 
the previous considerations. Accordingly, there must be an additional aspect of the 
geometry that has so far been neglected. Regarding Fig. [l](right) this issue is apparent: 
the pores are not straight but rather meandering. In consequence the length L of the 
path from any point A to another point B is always larger than the length D of the 
direct interconnection of the two points. To correct the pore length for the larger flow 
path an additional factor of ^ for the tortuosity must be introduced. Assuming r = 3.6 
this consideration yields for the Vycor pore network L ^ 1.2 D. This result can vividly 
be interpreted as follows: the shortest way from the bottom of the previously introduced 
sample cuboid to its top is about 20 % longer than its edge length a. 

2.1.2. Darcy's law With all the preceding considerations in mind one is able to derive 
an expression that describes the flow of a liquid through a porous network. For a given 
porous matrix with cross-sectional area A and thickness d (along which the pressure 
drop Ap is applied) the normalized volume flow rate ;^V^ is determined by 



-.- 



K 

rjd 



Ap . 



(4) 



This expression is also known as Darcy's law [14J. The proportionality constant K is 
the so-called hydraulic permeability of the matrix. It is given by 



K 



90 2 



(5) 



At this point it must be emphasized that the permeability is solely specified by the 
matrix' internal structure and consequently it should be independent of the liquid and 
of the temperature. 

2.2. Influence of the spatial confinement 

So far we have completely neglected that the mean pore diameters of the pore network 
investigated here are orders of magnitude smaller than characteristic in usual fiow paths 
in common miniaturized fiuid manipulating applications. Indeed, the pore radii are 
merely 10 to 100 times larger than typical molecular diameters of simple liquids like 





Figure 1: Illustration of the meaning of the tortuosity r of a pore network with 
isotropically oriented (left) and meandering (right) pores. 
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water. Within the systematic study of chain-hke hydrocarbons, which will be presented 

in section |5j the dimensions of the liquid's building blocks even exceed the channel's 

diameter. For that reason it is evident that some questions about the influence of the 

confinement on the fluid dynamics arise. In the following the two most apparent ones 

will be discussed. 

Though, beforehand, we will point out a remarkable feature being inherent in micro- 

and nanofiuidic devices and applications. Because of the tiny characteristic length scales 

C the Reynolds number 

inertial force pv C , , 

Re = 7 = (6) 

VISCOUS force rj 

(with the liquid's density p) usually fulfills Re <C 1. Thus, flows in such restricted 
geometries are most likely laminar rather than turbulent. From this point of view the 
implementation of the Hagen-Poiseuille law seems to be justified. But, for very high 



flow speeds v this assumption does not remain valid. In section 2.6 it will be shown that 



this fact plays a crucial role in the very initial phase of a capillary rise experiment. 

2.2.1. Validity of continuum mechanical theory Up to now we have assumed the law 
of Hagen-Poiseuille to be valid even in pores with diameters below 10 nm. However, 
one must not forget that this law is based on the principles of continuum mechanical 
theory, in which the behavior of a fluid is determined by collective properties such as 
the viscosity rj and the surface tension a. This assumption certainly holds for ensembles 
of 10^^ molecules. But within the pore confinement such amounts are not reached. This 
can easily be seen in the following example. Assuming water molecules to be spheres 
with a radius of 1.5 A in a hexagonal close-packed structure one arrives at only 1000 
molecules per cross-sectional area. As a consequence, the validity of the continuum 
theory has to be put into question. 

On this score especially the development of the surface force apparatus (SFA) 
has stimulated extensive studies over the last three decades. The mobility of water 
and several hydrocarbons in extremely confined films was examined by experiment 
[13 Uni [H] and in theory [IB]. These studies revealed a remarkable robustness of 
the liquids' fluidity down to nanometer and even subnanometer spatial confinement. 
Moreover the validity of macroscopic capillarity conceptions at the nanoscale was 
demonstrated [121 [20] • The measurements within this study will provide further hints 
whether the concepts of viscosity and surface tension still remain valid in nanopore 
confinement. 

2.2.2. Validity of the no-slip boundary condition The law of Hagen-Poiseuille implies 
the no-slip boundary condition. This means that the velocity of the fluid layers directly 
adjacent to the restricting walls equal the velocity of the walls themselves. Nowadays it 
is indisputable that this assumption does not hold unreservedly. Already 60 years ago 
Peter Debye and Robert Cleland introduced both slipping and sticking fluid layers at 
the pore walls in order to interpret a seminal experiment on liquid flow across porous 
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Vycor [H]. In that way, they were able to quantitatively account both for increased as 
well as for decreased measured flow rates (compared to the predicted ones) within their 
examinations of the flow of hydrocarbons through porous Vycor. 




Figure 2: Illustration of the possible boundary conditions along with the corresponding 
parabolic velocity profiles in a cylindrical tube with radius tq. Mass transport takes 
place only where the streaming velocity is different from zero, (left): The reduction of 
the net flow rate is due to sticking layers at the pore walls, which do not participate 
in the mass transport. In addition the maximum velocity in the pore center is smaller 
than for no-slip boundary conditions (middle) because of the smaller hydrodynamic pore 
radius r-^ < tq. This gives rise to a further dramatic decrease in the flow rate, (right): In 
contrast, a slipping liquid with a hydrodynamic pore radius rh > ro causes the highest 
streaming velocity and consequently the highest net flow rate. 

The concepts of a sticking and of a slipping liquid compared to the traditional 
no-slip boundary condition are exemplified in Fig. [2] for a cylindrical tube with radius 
tq. The degree of slip can be quantified by the slip length b with tq = rh — b. The 
hydrodynamic pore radius rh measures the distance from the pore center to the radius 
where the streaming velocity reaches zero. In this representation the sticking layer 
boundary condition is indicated by a negative slip length b whereas a positive slip 
length is typical of a slip boundary condition. The standard no-slip condition yields 
6 = meaning tq = r^. 

Up to date many factors have been found that seem to influence the boundary 
conditions. The most prominent and maybe the least controversially discussed 
amongst them is the fluid-wall interaction expressed in terms of the wettability 
[2Il|221|2ni[23[2Sl[2Sl[2Zll2Sll22]. The weaker the interaction is the more likely is 
slip. In addition, shear rates beyond a critical value are supposed to induce slip, too 
[301 ED E21 E3]- in contrast, the influence of surface roughness is rather debatable [M]. 
There are results for a decrease [551 122] as well as for an increase ^B] of the slip length 
with increasing surface roughness. Furthermore, dissolved gases [23 EH], the shape of 
the fluid molecules [20] or the add-on of surfactants [30] seem to influence the boundary 
conditions. To sum up, there is a huge set of factors (see Refs. [33 SOI E] for some good 
reviews) and certainly a complex interplay between many of them finally determines the 
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interfacial behavior. 

The permeabihty K of the membrane reflects such modified velocity boundary 
conditions: 

Srrg 8r rg ' ^' ^ 

Equation ([T]) illustrates the high sensitivity of K on b, provided b is on the order of or 
even larger than tq. Therefore, measuring the hydraulic permeability gives direct access 
to the slip length b for a given liquid under given conditions. 

One has to keep in mind that boundary conditions and fluid properties derived from 
measured flow rates are subject to a central restriction: one cannot verify the predefined 
parabolic shape of the velocity profile in the nanoscopic flow geometry. This is because 
there is no direct access to the profile itself but only to flow rates, which correspond 
to the velocity profile integrated over the whole pore cross-sectional area. Nevertheless, 
molecular dynamics simulations prove the formation of parabolic flow profiles even down 
to channel radii of three molecular diameters |121 SSI El] and, hence, justify inferences 
based on this major assumption. 

2.3. Laplace pressure & capillary rise 

From the physisist's point of view spontaneous imbibition is an impressive example for 
interfacial physics. The driving force behind the capillary rise process is the Laplace 
pressure pl acting on the curved meniscus of a liquid in a pore or porous structure. It 
is specified by 

with a being the surface tension of the liquid and -Ri and R2 being the local radii of 
curvature of the liquid surface. For that reason it is important to obtain information on 
the shape of the meniscus. In general it is determined by an interplay of surface forces 
and gravity. The influence of the latter can be estimated by means of the so-called 
capillary length 

Ac = ./^ (9) 

^ P9 

with g being the acceleration due to gravity and p representing the density of the liquid. 

For water (with a ~ 72 ^ and p = 1 ^^) Eq. to| yields a capillary length of 2.7mm. 

For surfaces with characteristic lengths smaller than Ac the gravitational force can be 

neglected. Therefore, in the porous structure of Vycor with mean pore diameters on the 

order of 10 nm the menisci are solely determined by surface forces. Consequently, they 

take the shape of a spherical cap (see Fig. [s]) described by a unitary radius of curvature 

i?i = i?2 = re resulting in a uniform Laplace pressure pi, = — at each point of the 

meniscus. 
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Figure 3: Sketches of cuts through pores before and during an imbibition experiment, 
(left): Spherical meniscus of a liquid in a pore with the static contact angle Oq. From 
simple geometric considerations one can deduce the relation between the radius of 
curvature r^ and the pore radius Tq: r^ = ^^^g . (middle): Water coating of the silica 
pore walls because of the finite humidity in the laboratory, (right): Illustration of the 
difference between the pore radius vq and the Laplace radius tl due to preadsorbed 
water layers, meaning that r^ 



COS 9q ' 



Eventually, from simple geometric considerations, one can deduce the relation bet- 
ween the radius of curvature Vc and the pore radius tq (see Fig. [3](left)) yielding 

2cr cos 6*0 



PL 



'"o 



(10) 



for the Laplace pressure. It is obvious that the interaction between liquid and matrix 
expressed in terms of the static contact angle 6o plays a crucial role in the capillary rise 
process. Spontaneous imbibition can only occur for a wetting or partially wetting liquid 
that is a liquid whose contact angle with the matrix is smaller than 90°. The influence 
of the contact angle hysteresis and a possible substitution of the dynamic contact angle 
6*0 for 6*0 will be discussed in section |2l6| 

Of course, as the liquid rises beyond its bulk reservoir to a certain level h the 
hydrostatic pressure pi^ = p g h acting on the liquid column increases. Accordingly the 
driving pressure Ap must be modified to Ap = p^ — Ph- The final state is then derived 
from a balance between pl and ph [IS]- The maximum rise level /imax is given by Jurin's 
law: 

2(7 cos^o 



hr 



11 



pgro 
It is worthwhile calculating /imax for water in the silica network of Vycor. Assuming 
6*0 = 0° for water on a glass substrate and vq ^ 5nm Eq. (11) yields /imax ~ 2.9 km 



corresponding to a Laplace pressure of 290 bar. Thus, for the here examined rise levels 
restricted by the maximum sample height to less than 5 cm the gravitational force can 
be neglected meaning that 

Ap = PL . (12) 
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Figure 4: Raytracing illustration of a Vycor sample during a capillary rise experiment 
filled up to the level h(t). 



The prevalence of surface forces over gravitation can likewise be expressed in terms of 

the dimensionless Bond number Bo: 

^ gravitational force PQ C? , , 

Bo = = ^ — < 1 13 

capillary force a 

because of the tiny characteristic length scales C. 

2.4- Influence of preads orbed liquid layers 

In order to gain information on the dynamics of the imbibition process we can now resort 
to Darcy's law (Eq. Q) discussed above. Considering Fig. |4] one is able to conclude that 
the sample height d appearing in Darcy's law has to be replaced by the actual rise level 
h{t) since only the filled parts of the sample contribute to the flow dynamics. Moreover, 
at any given time t the imbibed fluid volume V{t) is closely related to h{t) via 

V{t) = (j)iAh{t) . (14) 

Here 0i denotes the initial volume porosity of the Vycor sample. This porosity is reduced 
with respect to 0o because of water layers on the silica pore walls, which are immediately 
adsorbed under standard laboratory conditions because of the finite humidity and the 
highly attractive interaction between water and silica (see Fig. |3](middle)). This wall 
coating amounts to a 15 % to 20 % decrease in porosity. The adsorbed water is highly 
stabihzed and can only be removed at elevated temperatures |1B] . 

The exact degree of coating depends on the sample temperature and the absolute 
humidity, thus it cannot be known beforehand. Fortunately (pi can be extracted for each 
sample from the performed mass increase measurements, which will be presented later 
(see Fig. Isl). Given the density p of the imbibed liquid and the volume V^ of the sample 
block one gains direct access to the initial porosity via the overall mass increase M due 
to the liquid uptake (indicated in Fig. [s]): 

M , , 

0i = -rr • 15 

Similarly, the Laplace pressure must be influenced by the initial wall coating. The 
preadsorbed water layers necessarily lead to a reduction of the radius of curvature of the 
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menisci (see Fig. |3](right)). This effect will be taken into consideration by substituting 



Tl for Tq in Eq. (10) with Tl < Tq. We do not have any reliable information on the layer 
thickness x = tq — tl and hence on the Laplace radius tl. It can only be estimated from 
the difference between 0o and 0i assuming the pores to be perfect cylinders with a radius 
Tq. Then tl = tq -i/0i/0o enables one to estimate the upper bound x < 5 A (depending on 
temperature and humidity). Based on this result we will assume tl = (rg — 2.5 A) ±2.5 A 
in order to account for this effect. Fortunately, the impact of tl on the overall imbibition 



dynamics is comparatively small (see Eq. (20)) and, therefore, the uncertainty in the 
actual Laplace radius has only little effect on the final results. 



2.5. Dynamics of the imbibition process: BCLW law 

With all these preliminary considerations in mind we will now return to the dynamics 
of the imbibition process. Darcy's law Eq. ^ in conjunction with Eq. (14) leads to the 
simple differential equation 

K 

h{t) h{t) = -— Ap (16) 



(piT] 



solved through 



hit) 




ApVi . 



(17) 



Nowadays this square root of time behavior is often referred to as the Lucas- Washburn 
(LW) law after its alleged discoverers [lH HH]. Actually, the \/t-scaling was first found 
by Bell and Cameron ^\ more than 10 years earlier [50]. Hence we will term Eq. (17) 
BCLW law in order to acknowledge all contributors. 



Given the driving pressure Ap in accordance to Eq. ( 10 ) and the permeability K 



corresponding to Eq. ([T]) one obtains the rise level h at a certain time t 

a cos^o 



hit) 



2017/ 



r Vi . 



(18) 



Ch 




Figure 5: Schematic view of the capillary rise in a cylindrical pore with pore radius Tq, 
Laplace radius tl and hydrodynamic pore radius r^- 
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Along with Eq. (14) and the density p of the hquid Eq. (18) results in 



m(t) = pAf-l^rV-t, (19) 

^ ^ ^ 

Cm 

the mass increase m of the sample due to the liquid uptake as a function of the time 
t. The corresponding flow configuration is illustrated in Fig. |5j We introduced the 



imbibition strength T in Eqs. (18) and (19), which is given through 

r = i./^. (20) 

Comparable to the hydraulic permeability it depends solely on the matrix' internal 
structure. 



Interestingly, relation ( 20 ) reveals that in principle the capillary rise dynamics are 
proportional to the square root of the capillary's radius: F oc ^/r. As a consequence, 
in common capillary rise experiments with channel diameters on the order of some 
hundred micron the rise dynamics are so fast that without any additional instruments 
merely the steady-state equilibrium configuration at the very end of the process can be 
observed. This is contrary to imbibition in meso- and nanoporous networks. Here, the 
much smaller pores induce a lot more viscous drag, which reduces the overall dynamics 
significantly. Typical rise times in our experiments are on the order of some hours or 
even days. This fact allows for an easy recording of the dynamics by measuring the 
sample's mass increase m{t) as a function of the time. 

2.6. Remarks on short time effects 

Finally, we will outline some effects that dominate the imbibition process in its very 
initial phase, meaning the first few nanoseconds. However, these effects occur on time 
scales far beyond the temporal resolution limit of the experiment (on the order of one 
second) and therefore do not play a decisive role for the measurements. 

2.6.1. BCLW Law in Nanopores The BCLW law neglects inertial effects, thus it is only 
valid for time scales sufficient to establish viscous flow. For very short times a more 
general ansatz has to be applied (Bosanquet equation), which approaches the BCLW 
law for long times. The observed deviations from the BCLW law are possibly caused 
by a high degree of turbulences and initial deformations of the velocity profiles. This 
lasts until the fiow reaches its typical low Reynolds numbers, for which turbulent fiow 
cannot play a role at all. 

The typical time rinit that is required for the viscous fiow to establish itself in the 
pore can be estimated with j51j 



Tinit ^ ^ . (21) 

7] 

This typically yields times on the order of some 10 ps in high agreement with recent 
simulation studies 
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According to Ref. [HS] the fluid-wall interaction, and in particular the surface friction 
due to molecular corrugation, additionally seem to influence Tinit- For the atomic smooth 
CNT Tinit is rather increased and a lengthy non-BCLW behavior can be observed. The 
nearly frictionless surface of the nanotubes is also supposed to be responsible for the 
enhanced mobility of water and hydrocarbons flowing through these tiny ducts [S51 EZ] • 
Nevertheless, this phenomenon occurs on time scales far beyond the temporal resolution 
limit of the experiment and does not affect our measurements at all \5T[ \5E\ . 

2.6.2. Dynamic contact angle The static contact angle 9q of a drop resting on a surface 
is a constant determined by an interplay of the interactions between the solid, the liquid 
and the vapor phase. As a rule this assumption fails if the contact line begins to move. 
Here, a different contact angle, the so-called dynamic contact angle 9d is observed. 

It is this dynamic contact angle that is required as a boundary condition for 
modeling problems in capillary hydrodynamics. In consequence, it is important to 
estimate the changes that arise from this effect. Unlike the static contact angle, ^d is 
not a material property. Actually, for small spreading velocities v expressed in terms of 
the capillary number Ca 

Ca = viscous force ^ ^ r/ ^ 
capillary force a 

it seems to solely be influenced by the capillary number itself [SHI- This is the statement 
of the most widespread work relation describing the dynamic contact angle, namely the 
Hoffman- Voinov- Tanner law 6'^ oc Ca. It is valid for Ca < 10"'* [SH] or, with some 
correction, for Ca < 10~^ [20] ■ Interestingly, Richard Hofmann carried out a systematic 
study of dynamic contact angles in glass capillary tubes for a wide range of capillary 
numbers [^I] and some years later Jiang et al. gave an empirical correlation for his 
results [S2]: 

^"^^°:^"^^" = tanh (4.96 ■ Ca°-^°^) . (23) 

cos 6*0 + 1 ^ ^ 

This expression can be used to estimate the impact of the contact angle hysteresis 

on the capillary rise experiments presented here. For this purpose one needs information 

on the prevailing capillary numbers being tantamount to the advancement speed v of 



the liquid front. One can deduce this speed as the time derivative of Eq. (18): 



The obvious divergence for t — !■ s is again the manifestation of the deficiency of the 
BCLW law for very short time scales and has already been thematized in the literature 

[S31EI]. 

We obtain as the absolutely highest capillary number occuring in our experiments 



Camax ~ 10~^. Accordingly, the above approximation for low Ca Eq. (23) can be 
applied. With the static contact angle 9q = 0° this finally yields the highest achievable 
dynamic contact angle 6'D,max = 2.2° with cos6'D,max = 0.9993 compared to cos^o = 1- 
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Figure 6: Exploded view of the imbibition setup (raytracing illustration). The top 
of the Vycor sample is glued to a wire and hangs freely into the cell. The cell itself 
consists of a casing and a container. Both are build out of copper and can separately 
be thermostatted. 



It is apparent that the phenomenon of the contact angle hysteresis needs no further 
consideration. This result agrees with MD studies of Martic et al. |S1]. Indeed, they 
found an initial variation of ^d for liquids sucked into channels with 5 nm radius. 
Nevertheless, it relaxed within some 10 ns towards 6q what is again far beyond the 
temporal resolution limit of the experiment. 

What is more, the contact angle remaining 0° during the whole measurement entails 
an interesting side effect. In a recent lattice Boltzmann simulation Kusumaatmaja et al. 
[55] examined the influence of surface patterning on the capillary filling of microchannels. 
They found that for 6 < 30° the structuring has only little effect on the filling process 
and for completely wetting liquids there is no observable influence at all. Therefore we 
do not have to worry about atomic irregularities of the pore shape. 



3. Experimental 

Imbibition dynamics are studied via a recording of the samples' mass increase m{t). 
The corresponding setup is depicted in Fig. [6j For a time-dependent measurement of 
the force acting on and, hence, of the mass increase of the porous Vycor block the 
sample is installed on standard laboratory scales applying a special mounting. In order 
to perform measurements beyond room temperature a cell was constructed that allows 
for a simultaneous thermostatting of the sample itself and the liquid reservoir beneath. 
The samples were cut into regular shapes of either cylinders (V5) or cuboids (VIO) 
of known dimensions (cross-sectional area A and height ho). Prior to measurement they 
were stored in a desiccator to avoid the uptake of impurities from the surrounding air. 
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For highly volatile liquids such as water and the n-alkanes up to decane a sealing 
of the samples' top and side facets is of utmost importance. Partly this is because of 
the possible influx of molecules via the vapor phase. However, with increasing rise level 
another problem gains importance, namely the evaporation from the sample. This would 
actually lead to a pinning of the rise level h determined by a balance between the liquid 
supplying imbibition rate and the liquid detracting evaporation rate (see Fig. [?]). For 
water such a sealing can very well be realized with Scotch tape. For the short-length 
n-alkanes a two-component adhesive was applied. For all other liquids the samples' 
facets remained unsealed. 

The measuring principle of the gravimetric experiments can nicely be illustrated 
referring to the representative mass increase measurement depicted in Fig. |8| Four 
distinct regimes are indicated. In the beginning the sample is put into the thermostatting 
cell and is mounted on the laboratory scale so that it hangs freely above the bulk 
reservoir. Regime (a) is now described by a mass decrease of the sample, which can be 
seen as the obvious manifestation of its thermalization. Due to the elevated temperature 
(here T = 59 °C) the sample loses a part of its preadsorbed water and accordingly gets 
lighter. A sufficient thermalization is then given for an adequate convergence of the 
mass towards mg. 

The measurement is started by moving the cell upward until the sample touches the 
liquid surface. Immediately a meniscus forms along the block's perimeter. This induces 




Figure 7: Series of pictures of a pinning experiment of water invading a porous Vycor 
block (VIO, Vs = 8.14 x 6.05 x 29.00 mm'^) at room temperature. The sample is not 
sealed. Intervals between pictures: 106 min. Due to the increasing impact of evaporation 
the front is eventually pinned. 
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Figure 8: Representative measurement of the mass increase of a porous Vycor block 
(VIO, Vs = 6.25 X 5.05 x 12.70 mm^) due to the imbibition of n-tetracosane (n-C24H50) 
at T = 59 °C. Sohd hues are a y/i-Ht to the hquid's imbibition behavior and a hnear fit 
to the level of saturation. The arrow indicates the difference between the sample's mass 
before (mg) and after {me) the experiment. The experiment is divided into four regimes 
as indicated by the dashed lines and marked by (a) - (d). The data density is reduced 
by a factor of 2000. 



a force Fs towards the reservoir acting on the sample, which becomes noticeable by a 
sudden jump in mass right at the beginning of regime (b). The connection between F^ 
and the perimeter length P is determined by 

Fs = P(Jcos9o (25) 

widely known from the Wilhelmy plate assembly for measuring the surface tension of a 
liquid. Here it enables one to quantitatively analyze the mass jump. Assuming 6o = 0° 
and given the surface tension a ~ 30 ^^ of n-tetracosane at T = 59 °C Eq. (25) gives 
Fs ~ 0.7mN or, equivalently, an apparent mass jump of Am ^ 0.07g. This is in high 
accordance with the measurement. 

This jump is only an initial effect and can be seen as a constant offset that should 
not disturb the subsequent imbibition process itself. The latter is the outstanding effect 



in regime (b). It can be described by a \/t-fit in accordance to Eq. (19). Thus, the 



^/t-fit provides Cm and then Eq. (19) results in the imbibition strength 
p A V 0i o" cos Oq 



(26) 



where the initial porosity (/)[ according to Eq. (15) and the bulk fluid parameters are 
applied. The contact angle is assumed to be 0° and consequently cos 6*0 = 1. The latter 
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assumption was also tested experimentally and proved true for all liquids taken into 
consideration for this study by means of contact angle measurements. 

At some point a level of saturation is entered. The deviation from the -s/t-law 
in regime (c) signals that the sample is completely filled. In regime (d) the sample is 
detached from the bulk liquid reservoir in order to extract the mass rrie of the completely 
filled sample. The overall mass uptake is then M = rrie — rng. 

4. Spontaneous imbibition of water in porous Vycor 

Experiments on the capillary rise of water in V5 and VIO were conducted at three 
different temperatures T (25 °C, 40 °C, and 60 °C). The results are summarized in Tab. fl] 
As predicted for the imbibition strength the obtained values are independent of the 
temperature T. The arithmetic averages over all temperatures are (123.7±2.1) ■ 10~^ a/hi 



for V5 and (173.7 ± 3.1) ■ 10^'^ y^ ^^^ VIO. According to Eq. (20) these measured 
quantities correspond to hydrodynamic pore radii rh = (2.85 ± 0.21) nm for V5 and 
Th = (4.41 ± 0.28) nm for VIO, or, in terms of the slip length 6 = rt — rg, 

6 ~ — 5 A for water invading V5 and 
b^—5A for water invading VIO. 

Consistently for both V5 and VIO a negative slip length of about 5 A was found. 
Considering the diameter of a water molecule of approximately 2.5 A this result can 
be interpreted as follows: two layers of water directly adjacent to the pore walls are 
immobile meaning that they are pinned and do not take part in the flow. Taking this into 
account (through the above-mentioned negative slip length) the imbibition dynamics 
can conclusively be described by the model proposed in section |2j In consequence 
one may conclude that the residual inner compartment of the liquid obeys classical 
hydrodynamics based on continuum mechanical theory. As a result, our experiments 
confirm former findings on the conserved fluidity [151 dSl IE] and capillarity [ISl ED] of 
confined water - except for the pinned layers. 

Table 1: Characterization of the imbibition dynamics of water in porous Vycor at 
three different temperatures by means of the imbibition strength F as deduced from 



gravimetric measurements according to Eq. (26) (in units of 10 ^ v^)- 



temperature V5 VIO 

25 °C 126.7 ±8.2 170.0 ± 7.9 

40°C 120.5 ±7.4 180.9 ±11.6 

60°C 125.3 ±11.2 174.3 ±11.4 
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The assumption of the compartmentation of nanopore-confined water just stated 
is corroborated by recent molecular dynamics studies on the glassy structure of water 
boundary layers in Vycor and the expected existence of sticky boundary layers in Hagen- 
Poiseuille nanochannel flows for strong fluid- wall interactions [BSl EZl ESI ESI EDI EI] • By 
means of X-ray diffraction distortions of the hydrogen-bonded network of water near 
silica surfaces were found [72], which might be responsible for the markedly altered 
liquid properties. Tip-surface measurements document a sudden increase in the viscosity 
by orders of magnitude in 0.5 nm proximity to hydrophilic glass surfaces [75]. It also 
extends former experimental results with respect to the validity of the no-slip boundary 
condition for water /silica interfaces, in which this condition was proven down to at least 
10 nm from the surface ^1] whereas slip-flow for water is only expected at hydrophobic 
surfaces [73 El] • 

Here, it is important to highlight that a significant influence of silica dissolution in 
water can definitely be ruled out. The resultant aqueous electrolyte (often referred to 
as polywater) would of course be characterized by altered liquid properties and could 
therefore likewise be responsible for the observed decreased overall invasion dynamics. 
Nevertheless, as was pointed out in detail by Tombari et al. [TB] , the dissolution of silica 
in water for the measuring times and temperatures considered here is vanishingly small 
and may be neglected. 

Furthermore, any by-passing material transport via the vapor phase could be 
excluded. According to recently conducted lattice Boltzmann simulations such 
condensation processes would entail a breakdown of the usual imbibition dynamics 
[77]. However, a significant influence of evaporation from and subsequent condensation 
beyond the advancing liquid front can be ruled out referring to thermodynamic 
argumentations [TS]. Moreover, recently performed neutron radiography measurements 
did not only verify this prediction, but also allow us to conclude that the imbibition 
dynamics are only marginally affected by this transport mechanism [7U] . 

Finally, it is worthwhile noting that water encounters a negative pressure, which 
linearly decreases from — pl (on the order of some hundred bar) at the advancing menisci 
to atmospheric pressure at the bulk reservoir. Water's hydrogen bridge bond network 
is expected to be responsible for an increase in viscosity rj and decrease in density p 
under such large tensile pressures. Based on thermodynamic models for stretched water 
[50] I5T] , one may estimate a T-dependent ~ 3 % decrease in the overall dynamics due to 
this effect, which is, unfortunately, below the error margins. Finally, we would refer the 
interested reader to two more detailed publications reporting on capillary rise dynamics 
of water in nanoporous silica 



5. Spontaneous imbibition of hydrocarbons in porous Vycor 

The imbibition study of water in porous Vycor glass just presented reveals many aspects 
of the behavior of a liquid in interface-determined systems. Anyhow, these results could 
just as well have been obtained on a planar glass substrate, meaning that so far no 
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effects explicitly due to the liquids' restrictions to a cylindrical pore geometry have been 
reported. One might wonder whether the water molecules are too small as compared 
with the channel diameter to be subject to any limitations. This consideration motivates 
the hereinafter presented systematic study of the imbibition dynamics of chain-like 
hydrocarbons namely the homologous series of n-alkanes (n-C„H2„+2) and the branched 
hydrocarbon squalane. 

The alkane molecule's all-trans length £ is proportional to the number of C atoms 
n in the chain's backbone 

£(n) = 2.lA + (n-l) ■ 1.25A (27) 

whereas their diameter is independent of n (approximately 3 A). Hence, the alkane series 
allows one to easily vary the aspect ratio of the molecule, while many important liquid 
properties (density p, surface tension a, viscosity t]) nearly remain the same [HS]. In 
particular all alkanes totally wet silica. For that reason, alkanes are particularly well 
suitable for a systematic study of the influence of the shape of the liquids' building blocks 
on the flow dynamics through channels with diameters comparable to the molecule's 
length. 

Eleven linear hydrocarbons ranging from n-decane (n-CioH22) to n-hexacontane 
(n-C6oHi22) being tantamount to all-trans molecule lengths i between 1 nm and 
approximately 8nm were applied. Furthermore, the branched hydrocarbon squalane 
(2,6,10,15,19,23-hexamethyltetracosane) with six additional methyl side groups as 
compared to n-tetracosane (see Fig. [o]) was used. For better statistics the imbibition 
dynamics of each liquid/substrate combination were ascertained at least three times. 



The results of the study are summarized in Fig. 10 What is more, the invasion of 
n-hexadecane (n-Ci6H34) in V5 was also recorded as a function of the relative humidity 
(RH) of the surrounding atmosphere. For this purpose the whole setup was put into 
a box. Applying a humidifier the water content of the atmosphere within this housing 
could be increased to a maximum RH of approximately 50%. The resultant imbibition 
strengths F are listed in Tab.|2]together with the extracted values of the initial porosities. 
Obviously the invasion dynamics are suitably described through the arithmetic 




-v^*^ ^ 



Figure 9: Illustrations of the all-trans configurations of the linear hydrocarbon 
n-tetracosane (n-C24H5o) and the branched hydrocarbon squalane (C30H62). The 
backbones of both chains carry 24 C atoms, so they coincide in the all-trans molecule 
length ^ ^ 3 nm. 
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averages over all measured n-alkanes F = (122.8 ± 2.9) ■ 10^^ ^/m for V5 and F = 
(178.8 ± 2.5) ■ 10"''' ^/m for VIO. No distinct influence of the all-trans chain length i 
on the respective n-alkane's invasion dynamics can be detected. Only for the largest 
molecule n-hexacontane {i ~ 7.6 nm) in the smaller pores of V5 {tq = 3.4 nm) hints for 
significant deviations from the arithmetic average towards increased overall imbibition 
dynamics are found [F = (147.1 ± 5.0) ■ 10~^ V^]- Interestingly, enhanced dynamics for 
long-chain molecules (mostly polymers) have already been documented in the literature 
[m 1321 El] 3.nd could be ascribed to substantial slip at the liquid-substrate interface. 
This possibility will be addressed later on. 

Again, the hydrodynamic pore radii rh can be deduced from the measured 
imbibition strengths F. One attains rt = (2.85 ± 0.22) nm for V5 and rt = (4.49 ± 
0.28) nm for VIO, or in terms of the slip length b = r^ — Vq 

b ^ — 5 A for linear hydrocarbons invading V5 and 

b^ — 4A for linear hydrocarbons invading VIO. 

Comparable to the water measurements for both V5 and VIO a negative slip length of 
approximately 5 A is found. 

In this context it is worthwhile considering the impact of the pore walls' initial 
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Figure 10: Imbibition strengths F obtained from gravimetric measurements for a series 
of linear hydrocarbons (open symbols) ranging from n-decane up to n-hexacontane 
invading VIO (D) and V5 (O), respectively. Their mean values are visualized by the 
solid lines, lying within the confidence interval (dashed lines, confidence level: 90%). 
The corresponding all-trans molecule length £ of the alkane can be read from the top 
axis. The measured branched hydrocarbon squalane is indicated by the filled symbols. 
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Table 2: Imbibition strengths F and initial porosities (pi of n-hexadecane invading V5 as 
a function of the relative humidity (RH). 

RH (%) 0i r (10-^ y/m) 

24 0.275 116.5 ±3.8 

34 0.262 121.9 ±5.3 

50 0.245 121.9 ±4.1 



water coating on the liquid invasion. According to the results presented in Tab. |2] 
there is no significant influence of the RH on the overall rise dynamics, provided one 
takes into account the initial porosity as a function of the RH. Obviously the degree 
of the initial water coverage does not influence the hydrocarbon invasion. This result 
is especially important regarding Ref. [SS], in which already the presence of one to 
two water monolayers preadsorbed on a silica surface was found to be sufficient to 
considerably reduce the London force field of the glass. Nevertheless, the measured 
invasion dynamics prove that the modified alkane-substrate interactions do not change 
the wetting behavior of the hydrocarbon. 

However, especially the first adsorbed water layer is highly stabilized by the 
attractive potential between silica and water [IS]. In addition to that, a water coating 
on the pore walls is favored compared to an alkane coating. This is expressed by 
their different surface tensions and can be demonstrated by displacement measurements. 
Here, an initially full hydrocarbon filling in porous Vycor is replaced by a water filling 
through imbibition (after connecting the sample to a water reservoir). Cross checks of 
the sample's mass change with the different liquid densities reveal an absolutely complete 
substitution of water for the alkane. 

One has to assume that the imbibed hydrocarbon is not able to entirely displace 
the initial water coverage. At least the first layer adjacent to the pore walls should 
be present during the capillary rise process. This extremely immobile layer provides a 
first contribution of ~ 0.25 nm to the overall sticking layer thickness stated above. But 
there is still a residual of b, which can be ascribed to either a second water layer or 
to a layer of fiat lying hydrocarbon molecules. Thanks to their similar thicknesses one 
cannot distinguish between these two cases beforehand. 

Yet, the concept of a layer of fiat lying hydrocarbons is confirmed by the squalane 
measurements, which show decreased overall invasion dynamics. The ~ 15 % decrease 
in the imbibition strength and the consequent additional reduction of the hydrodynamic 
pore radii [rh = (2.61 ± 0.21) nm for V5 and r-^ = (4.19 ± 0.32) nm for VIO] by ~ 0.3 nm 
can uniformly be explained by the larger diameter of the squalane molecule as compared 
with a non-branched hydrocarbon chain. The extra methyl groups attached to the 
n-tetracosane backbone cause an increase of the sticking layer's thickness by a value 
comparable to the decrease in rh just mentioned. Interestingly, a similar infiuence of 
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the additional side branches on the boundary conditions was already reported [2S] • The 
continuously positive slip lengths b for both n-hexadecane and squalane revealed within 
this FRAP study notwithstanding, b was found to be always lower for the branched 
squalane. 

Summarizing, one may assume the sticking layer to be constituted of a water layer 
directly adjacent to the pore walls and a contiguous layer of flat lying hydrocarbon 
chains whereas the innermost liquid compartment obeys classical hydrodynamics based 
on continuum mechanical theory. Not until the chain length of the hydrocarbon reaches 
the diameter of the pore the results are seemingly influenced by the respective shape of 
the liquids' building blocks. 

The assumption of a sticky monolayer of hydrocarbons is confirmed by forced 
imbibition experiments on n-alkanes in Vycor [2], by studies regarding the thinning of 
n-alkane films in the surface force apparatus [H51 EH ESI ES] and by molecular dynamics 
simulations [SSIIHS]- On top of this. X-ray reflectivity studies indicate strongly adsorbed, 
flat lying monolayers of hydrocarbons on silica [^01 Ell 021 ES]- The retained viscosity 
of the confined alkanes was also verified by means of MD simulations |18j . 

With that the just stated results disagree with a series of recent publications. 
Especially the existence of a positive slip length b of short-chain alkanes was proven 
several times by experiment [221 [2S] as well as by simulation [Ml \TE\ , even for complete 
wetting and hence highly attractive liquid-substrate interaction. Interestingly, the 
hydrodynamic pore radius rh = (3.12 ± 0.26) nm of n-hexacontane in V5 reveals a slip 
length b = — 0.28nm, which is, admittedly, still negative. Nonetheless, it is markedly 
increased as compared with b = — 0.55nm for the arithmetic average over all n-alkanes 
by approximately 3 A - the thickness of the layer of flat lying hydrocarbons. This result 
might be interpreted as a step towards a slip boundary condition. Yet, because of the 
relatively large error bars it must be treated as a mere indication of altered boundary 
conditions. 

6. Spontaneous imbibition of silicon oils in porous Vycor 

So far the results on the flow of water and linear hydrocarbons in extreme spatial 
confinement revealed a retained fluidity and capillarity of the liquids. Considering the 
highly attractive interaction between fluid and substrate the verification of a sticking 
layer is not too surprising. Liquids slipping at a substrate are mainly observed in 
non- wetting configurations [751 HH 123 121] • Unfortunately such systems cannot be 
investigated by means of spontaneous imbibition, which always requires a wetting liquid. 
However, applying polymeric or oily liquids also seems to facilitate the occurrence 
of slip at the fluid-solid interface [22] • This presumption is corroborated by a recent 
imbibition study of polystyrene [SI]. What is more, the MD simulations carried out in 
Ref. [S] demonstrate a slip boundary condition for the capillary rise of a short-chain 
macromolecule (representative of the building blocks of silicon oil) with slip lengths on 
the order of 30 % of the tube radius. Therefore, when one tries to observe modified 
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Figure 11: Illustration of the molecular structure of l,l,3,5,5-Pentaplienyl-l,3,5- 
trimethyltrisiloxane. It is the main constituent of the applied Dow Corning diffusion 
pump oils DC704 and DC705. Five phenyl and three methyl groups are regularly 
attached to the trisiloxane (Si302-R8) backbone. The molecular dimensions can roughly 
be estimated using the diameter of a phenyl ring (~ 5.4 A) to be between Inm and 
1.5 nm. 



boundary conditions it will be rather promising to use such liquids. 

For this purpose we applied the diffusion pump oils DC704 and DC705 from Dow 



Corning. Their main constituent is the siloxane depicted in Fig. [TT] with molecular 
dimensions between Inm to 1.5 nm. They are characterized by a negligible vapor 
pressure and relatively high viscosities up to 200mPa-s. Because of the latter the 
capillary rise process is substantially slowed down. The BCLW-v^ behavior can be 
observed for weeks. Again, this preeminently elucidates the retained fluidity and 
Newtonian character of the liquid. The extracted imbibition strengths F are listed 
in Tab. [3l 

Compared to the values of water or the hydrocarbons the dynamics of the silicon 
oils are decreased. This is equivalently expressed in terms of the resultant slip lengths 
b (see Tab. [s]). They all indicate a sticking layer boundary condition whose thickness is 
approximately 1.4 nm. Interestingly this result is independent of both the sample type 
and the silicon oil. Thus it is not too erroneous to conclude that, again, one layer of 

Table 3: Characterization of the imbibition dynamics of Dow Corning silicon oils in 
porous Vycor at room temperature by means of the imbibition strength F (in 10~^ ^/^) 
and the resultant slip length b (in nm). 

DC... X^ Xl° 

T b T b 

704 63.6 ±5.6 -1.35 ±0.23 109.2 ± 12.6 -1.39 ±0.40 

705 53.2 ±4.8 -1.52 ±0.21 122.4 ±11.4 -1.19 ±0.38 
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siloxane molecules is pinned to the pore walls and does thereby substantially lower the 
invasion dynamics. 

This result does not necessarily have to contradict the findings of the studies 
mentioned at the beginning jUJ [321 El]- The applied silicon oils presumably are not 
chain-like enough but rather block-like. Possibly this internal molecular structure is 
decisive for the occurrence of slip at the interface. In this regard one has to recall the 
increased dynamics of n-hexacontane in V5 compared to the short-chain hydrocarbons. 
This result stimulates future studies on more extended alkane and polymer melts with 
molecular weights beyond 1 ^. 

7. Phase transitions of liquids in nanopores explored by spontaneous 
imbibition 

In this final section the infiuence of the nanopore confinement on temperature induced 
phase transitions will be examined. Usually phase transitions are accompanied by unique 
variations of the fiuid properties. Since the overall rise dynamics sensitively depend 
on these quantities one can easily detect such characteristic deviations. According to 



Eq. (19) the desired information is comprised in the prefactor Cm of the observed \/i 



behavior of the measured mass increase: 



C„^,/|.^^*f^r. (28) 

with only the quantities in the first factor depending on the temperature. Here, the 
slight T dependence of 0i can be neglected in the limited T regions observed. Moreover, 
the static contact angle 6q is zero for all T. The matrix properties are T-invariant 
and can be eliminated through normalization. Hence, by performing several imbibition 
measurements at different temperatures in the T-range of interest a simple comparison 
of the prefactors might already provide an insight into the phase transition behavior of 
the confined liquid. 

7.1. Meniscus freezing in n-tetracosane 

At first we present a T-dependent study of the capillary rise of the linear hydrocarbon 
n-tetracosane (n-C24H5o) in porous Vycor. The chosen alkane exhibits a surface freezing 
transition at Ts = 54 °C. Surface freezing is the formation of a single solid monolayer 
fioating on the bulk melt in a T-range between the bulk freezing temperature Tf and a 
temperature Tg [23 ESI EZl EHl EH] • Upon surface freezing the hydrocarbon molecules are 
rectified, parallel-aligned with their long axis along the surface normal, and the center 
of mass lattice is hexagonal, resulting in a 3 nm thick surfactant-like monolayer. The 
abrupt onset of surface freezing at Tg allows one to switch on (and off) the 2D crystalline 
phase at the free surface by a small T variation. The goal of this study is to investigate 
whether this particular molecular rearrangement is detectable by and how it affects the 
imbibition dynamics of the alkane in nanopores. 
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The surface freezing transition is accompanied by a distinctive signature in the T 
behavior of the surface tension a, known from bulk n-alkanes, that is a change from 
a small negative, above Tg, to a large positive T-slope, below Tg (see lower panel in 



Fig. 12) |96j. This distinctive behavior of the surface tension must be mapped in the 



T-dependence of the measured mass uptake. 

For convenience the measured prefactors Cm{T) of the i/t-imbibition dynamics 
are normalized by the prefactor value at an arbitrarily chosen temperature T'^. This 
procedure finally yields the normalized imbibition speed v^ of the measurement at the 
temperature T. One should not be confused by the term 'speed', since f" is actually 
dimensionless. Nevertheless, it quantifies the mass uptake rate of the sample at a given 




temperature T (°C) 



Figure 12: (upper panel): Measured normalized imbibition speeds v"^ (for T° = 60 °C) 
of tetracosane in V5 (D) and VIO ( _ ), respectively, in comparison to values calculated 
under the assumption of meniscus freezing (— ). Insets: Illustration of a surface frozen 
layer at an advancing meniscus (T < Tg) and of a sticky, flat lying boundary layer at the 
silica pore wall in the entire T- range investigated, (lower panel): T-dependent surface 
tension of tetracosane. Insets: Illustrations of surface freezing at a planar tetracosane 
surface. The dashed lines correspond to extrapolations of the calculated imbibition 
speeds and of the surface tension below Tg in the absence of surface freezing, respectively. 
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Figure 13: Illustration of the rod-like molecular structure of the liquid crystal 



octyloxycyanobiphenyl (80CB) with the phase transition behavior: crystalline — -V A 
— y N — y I. The molecule's head (comprising the cyano and the biphenyl group) and 
its tail (the octyl chain) each have a length of slightly more than 1 nm. Its overall length 
is ~ 2 nm. As opposed to the hydrocarbons the 80CB molecule is rather rigid. 



temperature. According to Eq. (28) this quantity can also be calculated: 



v^{T) 



C^{T) _ p{T) a{T) ■ r^{T-) 



(29) 



It is important to note that in Eq. (29) the geometry as well as the internal topology 



of the substrate do not play a role any longer. Assuming the (j{T) kink-behavior 



according to Fig. 12 along with the T dependency of r] and p one can calculate theoretical 
values of v^{T). They are plotted as a solid line in the upper panel of Fig. 12 The 



measured speeds are indicated by the single points. Due to the remarkable coincidence 
of the temperature characteristic of the kink in v^ and the degree of deviation from a 
suppressed phase transition we feel encouraged to solely attribute the distinct change in 
the imbibition dynamics to a change in cr(T) at the advancing menisci typical of surface 
freezing. The corresponding nanoscopic flow configuration is illustrated in the insets of 
Fig. 12 for both T < Ts and T > T^. We would refer the interested reader to a more 
detailed publication reporting on the surface freezing of n-tetracosane in nanoporous 

silica [inn]. 



7.2. Flow of a rod-like liquid crystal (80 CB) in Vycor 

Finally, the phase transition behavior of a thermotropic liquid crystal in nanopore 
confinement was examined. For this purpose we applied one of the best scientifically 



studied liquid crystals, the rod-like octyloxycyanobiphenyl (80CB, see Fig. 13 for an 
illustration of its molecular structure). At ~ 80 °C it undergoes a transition from the 
isotropic to a nematic phase, which is accompanied by the occurrence of a characteristic 
shear viscosity minimum. Again, this distinctive behavior must affect the rise dynamics 
which enables one to examine the phase transition in nanopore confinement. The most 



common phases are illustrated in Fig. 14 
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Figure 14: Illustration of common liquid crystal phases. Upon heating a frozen liquid 
crystal one might pass the following mesophases: smectic C (C), smectic A (A), and 
nematic (N). Beyond the clearing point T^ the liquid crystal is in the isotropic state (I). 



7.2.1. Nemato-Hydrodynamics: Shear viscosity minimum and presmectic divergence of 
flowing nematic liquid crystals The momentum transport in nematic liquid crystals 
shows an anisotropy since it depends on the mutual orientations of the macroscopic 
molecular alignment (the director n), the flow velocity {v) and the velocity gradient 
(V v). In 1946 Miesowicz defined three principal shear viscosity coefficients of nematics 
[TOlj . which can be measured in three different Couette flow experiments sketched in 



Fig. 15 Typically magnetic fields are applied in order to align the molecules in the 
nematic sample. Intuition suggests that the lowest resistance to the nematic fiow, i.e. 
the lowest viscosity value, should be 772- Among the two remaining viscosities, rji should 
have the highest value. 

Indeed, when the orienting magnetic field, i.e. the director n, is parallel to the 
velocity v of the nematic fiow, the lowest viscosity is recorded (see open symbols in 



Fig. 16). Nevertheless, this relatively simple picture of the viscosity of nematic liquid 



crystals is disturbed for the compounds exhibiting the transition to the smectic A phase. 
Then, with decreasing temperature, the viscosity 772 shows a strong increase and goes 
to infinity at the N-A phase transition. The rji and 773 viscosities are almost unaffected. 
At a temperature that is a few degrees below the temperature at which the transition 
to the smectic A phase takes place, the viscosities 772 and 773 interchange their roles and 
then the lowest nematic viscosity corresponds to the fiow in the 773 configuration. The 
presmectic behavior of the 7/2 viscosity is due to the formation of precursors of smectic 
planes with z\\v that would be immediately destroyed by the velocity gradient, thus 
this configuration is rendered unfavorable. 

The behavior of the freely fiowing compound obeys a general principle that can be 
formulated as follows: a free fiuid adopts such a manner of fiow, as corresponds to the 
minimum of its viscosity at given conditions |1U3] . Accordingly, the transition from the 
isotropic to the nematic phase manifests itself in a strong decrease of the shear viscosity 
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Vi 







Figure 15: The experimental Couette flow conditions for measurements of the three 
Miesowicz shear viscosity coefficients of nematic liquid crystals: r]i for n J^v and rT 1 1 V f , 
?72 for n\\v and n ± V f , 773 for n ± i7 and n _L V f . 
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Figure 16: Miesowicz shear viscosities 772 (□) and rj^ (A) of the liquid crystal 80CB 
compared to its free flow viscosity ?7iso ( • ) according to |lU2j . 



?7iso that is very close to r/2. Consistently, beyond the presmectic cross-over of ri2 and 
rjs the viscosity of the freely flowing liquid crystal ?7iso follows 773 (see ffiled symbols 



in Fig. 16). This result is interpreted in terms of rearrangements of the molecular 



alignment n with respect to the velocity field v, which can easily be assessed by means 
of examinations of the compound's viscosity. 

This distinctive behavior of the viscosity must be mapped in the T-dependence of 



the measured mass uptake. In Fig. 17 some representative measurements of 80CB 
invading V5 are shown along with their corresponding Vt-fits. In analogy to the 
procedure employed for the evaluation of the surface freezing transition we conducted a 
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quantitative analysis of the phase transition behavior of 80CB in nanopore confinement 
referring to the prefactors obtained from the fits. 

Again, the prefactors are normahzed by the prefactor value at an arbitrarily chosen 
temperature T^. The thereby obtained normalized imbibition speeds u'^(T) calculated 
from the fitting parameters are indicated by the single points in the upper panel of 



Fig. 18 Assuming the riiso{T) behavior according to Fig. 16 along with the T dependency 



of 0" and p one can calculate the theoretical values of ^"(T) based on Eq. (29). They 
are plotted as a solid line in the upper panel of Fig. [Tsj 

The results show a variety of astonishing features. First of all, both V5 and 
VIO reveal comparable characteristics of the invasion dynamics of the liquid crystal 
80CB. Anyhow, only in the nematic phase the T-dependent behavior of the measured 
imbibition speeds coincides with the prediction based on the ?7iso values presented in 



Fig. 16 In particular, the distinctive bump in the proximity of the clearing point Tc as 
the direct manifestation of the shear viscosity minimum in the theoretical behavior of 
v^ is unambiguously absent. The imbibition speed rather increases monotonously with 
increasing temperature. 

As mentioned before this distinctive feature of the viscosity at the clearing point 
is caused by the inset of an alignment of the molecules (in the nematic phase) with 
respect to the flow direction, that is n 1 1 -y. Therefore its absence can intuitively be 
interpreted in terms of an already existent alignment of the molecules beyond Tc. It 
is obvious to conclude that such an alignment is easily induced by the extreme spatial 
confinement to cylindrical pores with diameters that are not more than five times the 
length of the molecule itself. In addition, considering their rigidity and, in particular. 
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Figure 17: Specific mass uptake of V5 due to the imbibition of the hquid crystal 80CB 
as a function of the time for selected temperatures below and above the clearing point 
Tc ~ 80 °C. Solid lines correspond to v^-fits. The data density is reduced by a factor of 
2500. 
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Figure 18: (upper panel): Measured normalized imbibition speeds w° (for T^ = 75 °C) of 
80CB in V5 (A) and VIO (D), respectively, in comparison with values calculated on the 
basis of the viscosity values in the lower panel (— ). (lower panel): T-dependent viscosity 
of 80CB by way of comparison (see Fig. [16] for a detailed view). The dashed lines 
correspond to extrapolations of the calculated imbibition speeds and of the viscosity in 
the absence of the shear viscosity minimum and the presmectic divergence, respectively. 



the parabolic flow profile established in the pore, it is hard to think of any alternative 
to the tendency of molecular alignment parallel to the pore axis and, consequently, to 
the flow velocity v. From this point of view the absence of the shear viscosity minimum 
and, consistently, of the nematic to isotropic phase transition is not surprising at all 
but rather consequent. In this context it is more suitable to label the phase beyond the 
clearing point not isotropic but paranematic (P). This term preeminently discloses the 
solely confinement-induced alignment of the liquid crystal. 

It was demonstrated experimentally |1U4[ 1105^ I1U6L I1U7L 1108] . in agreement with 
expectations from theory |1U9[I11U]I111] . that there is no 'true' I-N transition for liquid 
crystals confined in geometries spatially restricted in at least one direction to a few 
nanometers. The anchoring at the confining walls, quantified by a surface field, imposes 
a partially orient at ional, that is, a partially nematic alignment of the confined liquid 
crystal, even at temperatures T far above the clearing point Tc. The symmetry breaking 
does not occur spontaneously, as characteristic of a genuine phase transition, but is 
enforced over relevant distances by the interaction with the walls. 

The assumption of such a paranemtic phase is corroborated by recent birefringence 
measurements of liquid crystals confined to an array of parallel, nontortuous channels of 
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10 nm mean diameter and 300 pm length in a monolithic silica membrane |112] . These 
measurements elucidate that the surface anchoring fields render the bulk discontinuous 
I-N transition to a continuous P-N transition. The transition temperature Tc is found to 
be changed only marginally, due to a balance of its molecular alignment induced upward 
and its quenched disorder (attributable to wall irregularities) induced downward shift. 
This agrees with the observations of liquid crystals imbibed in tortuous pore networks 
[TMfTIIS] . 

Interestingly, due to the complete absence of the shear viscosity minimum in 
the results of the T-dependent series of imbibition measurements shown in Fig. 



definition of a P-N transition temperature in confinement is not possible at all. This 
is confirmed by a simple extrapolation of the 772 viscosity to higher T, which saliently 
reproduces the measured imbibition speeds beyond the bulk clearing point. This can 
only be interpreted in terms of an extremely high degree of orientation already existent in 
the paranematic phase; at least higher than suggested by the birefringence measurements 
|112] . This difference is most probably due to the basically differing detection method 
of the P-N transition: normally birefringence or calorimetry (DSC) measurements 
performed with the confined static liquid are applied for this purpose. But, the viscosity 
measurements presented here refer to the liquid's dynamics in nanopore confinement. 
As already mentioned before, the additional emerging fiow velocity and in particular 
the velocity gradient seemingly enhance the paranematic orientational alignment in the 
772- configuration significantly. 



The second remarkable feature of the results presented in Fig. 18 is the absence of 
the presmectic divergence of the viscosity of the freely fiowing liquid crystal. This would 
result in a dramatic drop of the imbibition speed due to the inset of smectic layering. 
However, the measured values do not indicate such an effect. Its absence rather suggests 
a suppression of the A phase in favor of the N phase. This is elucidated by a simple 



extrapolation of the viscosity to lower T as indicated in Fig. 18, which preeminently 
reproduces the measured imbibition speeds. 

What are the reasons for this discovery? First of all, the confinement to a cylindrical 



pore (rather than to a film geometry like in the Couette fiow depicted in Fig. 15) renders 
the 7]3 viscosity as unfavorable as the r/i viscosity. This again clarifies the high stability 
of the nematic 7^2- configuration in confinement. A cross-over behavior as ascertained for 
the freely fiowing bulk liquid can hence be excluded. Yet, the presmectic divergence of 
ri2 is caused by the destruction of precursors of smectic planes for n\\v. From this point 
of view the suppression of the A phase is not surprising at all but a mere consequence 
of the overall stabilization of the r72-configuration in the nanopore confinement. 

Even the static nanopore-confined liquid crystal shows such modified mesophase 
behavior |106] . For example the heat capacity anomaly typical of the second-order N-A 
transition in rod-like liquid crystals immersed in aerogels is absent or greatly broadened 
|113t I114J . Nuclear magnetic resonance (NMR) measurements revealed the lack of 
pretransitional smectic layering due to the rough surface of the confining walls |115] . 
Furthermore, a systematic study of the infiuence of the degree of confinement indicates 
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that the N-A transition becomes progressively suppressed with decreasing pore radius 
whereas the stabihty range of the nematic phase is increased [116J . 

Finally, we will present some results from a more quantitative analysis of the 
measurements. According to the previously used evaluation method one arrives at slip 
lengths of 6 = (—1.11 ± 0.23) nm for V5 and b = (—1.54 ± 0.31) nm for VIO. Again, a 
sticking layer boundary condition has to be applied. 

In conclusion the T-dependent imbibition measurements of the liquid crystal 80CB 
revealed that confinement plays a similar role as an external magnetic field for a 
spin system: the strong first-order I-N transition is replaced by a weak first-order 
or continuous paranematic to nematic transition, depending on the strength of the 
surface orientational field |117] . Based on detailed knowledge of the static (equilibrated) 
liquid's behavior in the nanopores as deduced from previously conducted birefringence 
experiments, we were able to procure complementary results with respect to its dynamic 
(non-equilibrium) behavior. The additional emerging flow velocity and in particular 
the velocity gradient enhance the paranematic orientational alignment significantly 
rendering the P-N transition even broader than known from the equilibrium state. Due 
to the high stabilization of the 772- configuration in the N phase the A phase is suppressed 
and the stability range of the nematic phase is increased. Nevertheless, despite the 
molecular alignment no indications of velocity slip were found. 

8. Conclusions 

We have presented a gravimetric study on spontaneous imbibition of a series of molecular 
liquids as a function of the complexity of their basic building blocks. By thoroughly 
accounting for the influence of preadsorbed liquid layers this method allowed us to 
gain insights regarding the flow properties and the phase transition behavior of these 
systems in extreme spatial confinement. We inferred from our studies an interfacial 
boundary layer adjacent to the pore walls with a defined thickness whose dynamics are 
mainly determined by the interaction between liquid and substrate. This manifests 
itself in terms of a negative velocity slip length for all but a few (more polymer- 
like) molecular systems investigated here. The flow properties of the pore-condensed 
molecules in the pore center are, however, remarkable robust, that is bulk-like. This 
finding of two distinct species (with regard to flow dynamics) is reminiscent of the 
partitioning of pore condensates found in vapour sorption isotherms (film-condensed 
versus capillary-condensed state) |118] . of the two species often reported in studies on 
self-diffusions dynamics (slowed-down or increased self-diffusion versus bulk dynamics) 
[119^ 1120^ 1121^ I122J . of the effects of pore condensates on the deformation of the 
matrix (expansion versus contraction of the pore walls upon change from film-condensed 
to capillary-condensed state) |123] . Last but not least, it reminds of the distinct 
crystallization behaviour of pore condensates (an amorphous boundary layer versus a 
crystallized fraction of molecules in the pore center) |118j . In principle, this partitioning, 
but also the highly stabilized nematic phase of a confined liquid crystal, can be traced to 
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the strong influence of tfie substrate potential on the liquid or solid layers, respectively, 
right adjacent to the pore walls [124^ 112511126] . Thus, from a more general perspective, 
our study highlights that this interaction is not only of importance for an understanding 
of the equilibrium properties of nanopore-confined condensed matter |127lll28j . but also 
for the non-equilibrium, flow properties of the molecular assemblies studied here. 

For the future the investigation of even more complex molecules and macromolecules 
such as proteins in mesopores may be envisioned and of interest [129^ ^^ . Moreover a 
variation of the substrate species, e.g. a change from porous silica to readily available 
porous silicon |13U] . porous alumina, porous gold |13H I132[ 1133] or carbon nanotube 
[STf Eni 1134] membranes, and a modification of the pore wall chemistry, and the 
resulting change in the liquid/substrate interaction could lead to quite new imbibition 
phenomenologies. For such studies the concepts and findings presented here may set 
the stage for a better understanding. 
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